Study of two-neutron radioactivity in the decay of 2e O 
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A new technique was developed to measure the lifetimes of neutron unbound nuclei in the pi- 
cosecond range. The decay of 26 0— > 24 + 71 + 71 was examined as it had been predicted to have an 
appreciable lifetime due to the unique structure of the neutron-rich oxygen isotopes. The half-life 
of 2e O was extracted as 4.5ti's (stat.) ± 3(sys.) ps. This corresponds to 26 having a finite lifetime 
at an 82% confidence level and, thus, suggests the possibility of two-neutron radioactivity. 

PACS numbers: 21.10.Tg, 23.90+w, 25.60.-t, 29.30.Hs 
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Drastic changes in the structure, properties, and avail- 
able decay-modes of isotopes with extreme neutron-to- 
proton ratios have been observed in comparison to their 
stable counterparts [Irlij]. In I960 Goldansky predicted 
that the unique properties of very proton-rich nuclei 
would produce scenarios in which one- and two-proton 
radioactivity could be observed Q. These exotic modes 
of radioactivity were later verified through measurements 
of the one-proton decay of 151 Lu [jl and 147 Tm Q and 
the two- proton decay of 45 Fe [ill Till]. In particular, a 
measurement by Miernik et al. allowed for the exotic 
two-proton radioactivity of 45 Fe to be fully characterized 
for the first time providing insight into the three-body 
structure of this proton dripline nucleus fl2j . New modes 
of radioactivity have also been observed for neutron-rich 
nuclei including B- delayed two- [l3| , three- [3] , and four- 
neutron decays [15| . 

Recently, Grigorenko et al. calculated the lifetimes 
of one-, two-, and four-neutron decays from unbound 
nuclei While long lifetimes for proton decay can 

arise from the presence of both the Coulomb and an- 
gular momentum barriers, the neutron decay must only 
overcome the angular momentum barrier and will there- 
fore have a much shorter lifetime. However, Grigorenko 
et al. showed that it would be possible for some neutron- 
rich unbound nuclei to have long lifetimes, reaching the 
limit of radioactivity, depending on the nuclear struc- 
ture [l6|. While an exact limit for a lifetime to be con- 
sidered radioactivity does not exist, different arguments 
have been presented suggesting a lower limit between 
10~ 14 s and 10~ 12 s (see discussion in Refs. Q and 
For the two- neutron unbound 26 O a lifetime on the order 
of 10~ 12 s was predicted for a ground-state resonance en- 



ergy of about 150 keV and a pure v(d^/2) 2 configuration 
for the valence neutrons pj| . 

The recent measurement of the 26 O ground state res- 
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keV by Lunderberg et al. 
opened up the exciting possibility for a new type of ra- 
dioactivity to be discovered 17j. This limit was con- 
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firmed by the GSI-LAND group which determined the 
26 O ground state to be unbound by less than 120 keV [l8|. 
Upper-limits have already been placed on the lifetime of 
26 O from previous experiments. Searches for 26 O using 
fragment separators placed a lifetime limit of < 200 ns 
(roughly the flight time through the separator) based on 
the non-observation of the bound nucleus [l|| [2(| ■ The 
GSI-LAND group improved that limit by an order of 
magnitude to < 5.7 ns based on the reconstruction of 
an 26 fragment which decayed in flight outside the tar- 
get [IH- The 26 lifetime is therefore too short to be ob- 
served by traditional implantation/decay methods, where 
lifetimes down to several hundred nanoseconds have been 
achieved using digital electronics [2ll ] , and must be mea- 
sured in-flight. For two-proton emission, reconstruction 
of the decay vertex by Mukha et al. [22j and an adaption 
of the Recoil Distance Doppler Shift (RDDS) technique 
by Voss et al. [23| were used to measure the lifetime in- 
flight of 19 Mg. 

In this letter, the lifetime of 26 is extracted using a 
novel technique based on the velocity difference of the 
emitted neutrons and residual charged fragment travel- 
ing through a thick target. This approach is a variance on 
the Doppler-Shift Attenuation Method (DSAM) in which 
an excited nucleus is slowed in a solid material resulting 
in a distribution of Doppler shifted '7-ray energies that 
can be related to the lifetime [2J-[2g]. This concept is 
extended in the present work for nuclei which decay by 
neutron emission. The analysis is completed using the ex- 
perimental data from the work of Lunderberg et al. [l7| , 
in which the 26 O ground state resonance was originally 
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FIG. 1. (Color online) Experimental 24 + n + n decay energy 
spectrum (solid black points) with causality cuts applied is 
compared with the Monte Carlo simulation (solid black line) 
with two components: (red long-dashed line) the 26 ground 
state resonance and (green short-dashed line) the first excited 
state. The vertical dotted line represents the selection of 2e O 
events used in the analysis. 



measured. 

Since the experimental details have been provided pre- 
viously in Ref. [l7| , only a brief overview is presented. 
The 26 O was produced using a one-proton knockout re- 
action from a 82 McV/u 27 F beam produced at the Na- 
tional Superconducting Cyclotron Laboratory at Michi- 
gan State University. To produce the 27 F beam, a 
140 MeV/u primary beam of 48 Ca bombarded a 1316 
mg/cm 2 9 Be production target. The A1900 fragment 
separator [27} was used to select the desired 27 F frag- 
ments which were then impinged on a 705 mg/cm 2 9 Be 
reaction target in the experimental vault. The 26 0^ 24 
O+n+n decay was measured using the Modular Neutron 
Array (MoNA) and the 4 Tm superconducting dipolc 
magnet [28|. The dipole magnet bent the charged frag- 
ments about 43° into a suite of charged particle detectors, 
which allowed for the mass, charge, kinetic energy, and 
angle of the charged particle to be reconstructed from its 
track through the magnet [29|. MoNA was placed 6.05 m 
from the reaction target and provided the measurement 
of the velocity and angle of the neutrons. 

The three-body decay energy of the 24 + n + n system 
was calculated as .Edecay = M26Q — M24 — 2M„, where 
A/26 (M 2 4 ) is the mass of 26 ( 24 0) and M n is the 
neutron mass. The invariant mass, M260, was calculated 
from the experimentally measured four-momenta of the 
24 and two neutrons. The three-body decay spectra re- 
quires a triple coincidence of two interactions in MoNA 
that pass the causality cuts and a 24 O fragment. The 
causality cuts are used to select true 2n events from mul- 
tiple scatteringof a single neutron and are discussed in 
detail in Ref. [17| . 

A detailed Monte Carlo simulation was used to fit the 
experimental spectrum as described in Ref. (l7j . The 
simulation included all relevant components of the exper- 
imental setup. In particular, special care was taken in re- 
producing the neutron interaction observables in MoNA 
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FIG. 2. The decay of 26 within the thick 9 Be target is illus- 
trated for two cases: (top) very short lifetime corresponding 
to an immediate decay and (bottom) a lifetime around 30 ps 
which allows the 26 O to exit the target before decaying. 



using the Geant4 framework with the custom neutron 
interaction model MENATE_R (30[- As shown in Fig.Q] the 
experimental three-body decay spectrum was well repro- 
duced by the Monte Carlo simulation including the decay 
from both the ground state (red long-dashed line) and 
first excited state (green dashed line). It is important 
to note that the ground state resonance was determined 
from a fit of the data jl?} while the placement of the 
first excited state was taken from predictions from the 
continuum shell- model [31 1. 

Two scenarios for the decay of 26 O with different life- 
times are illustrated in Fig. [2J The 27 F beam enters the 
3815 /im (705 mg/cm 2 ) 9 Be target followed by the one- 
proton knockout reaction producing the 26 O. If the reac- 
tion was to occur at the beginning of the target where 
the 27 F beam is traveling at about 11.8 cm/ns and the 
26 had a very short lifetime (top of Fig. [2]) then the 
neutrons would be emitted with an average velocity of 
11.8 cm/ns. In the other case if the 26 had a lifetime 
of 30 ps (bottom of Fig. [2]), which is roughly the time 
of flight through the target, then the neutrons would be 
emitted with an average velocity of 10.9 cm/ns due to 
the energy loss of the 26 O fragment traveling through the 
target. Thus, the observation of a shift in the expected 
neutron velocity can provide a measure of the lifetime of 
26 0. 

The relative velocity between the neutrons and frag- 



ment is defined as V re i — V n — Vf rag , where V n (Vfrag) is 
the velocity of the neutron (fragment) in the laboratory 
frame. The relative velocity was examined to remove 
the effect of the momentum dispersion of the 27 F beam 
(Ap/p = 2%). Thus, the variation in the incoming ve- 
locity of the 27 F is removed event by event. Since the 
reaction point in the target is unknown, the fragment ve- 
locity {Vfrag) is calculated assuming the reaction occurs 
at the center of the target. The width of the V re i dis- 
tribution will be dependent on the target thickness and 
magnitude of the decay energy (both of which will in- 
crease the width). If the reaction point was known on an 
event-by-event basis and the decay energy was very small 
then the V re i distribution should be narrow and centered 
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around zero. While the width of the experimental V re i 
distribution will be increased, the mean velocity will still 
be centered around zero if the 26 O lifetime is short. Thus, 
a shift in the relative velocity away from zero would in- 
dicate a long-lived component of the decay. 

The 26 O events were selected from events passing the 
causality cut criteria and having a Edecay < 1.0 McV, as 
indicated from the dotted grey line in Fig. [T] This se- 
lection should maximize the statistics and minimize the 
contamination of other decay channels. Based on the fit 
of Fig. [T] the 26 ground state resonance accounts for 
96% of the events with E decay < 1.0 McV. The rela- 
tive velocity between the 24 and each of the emitted 
neutrons is shown in Fig. Eta). The experimental V re i 
(solid black points) is shifted away from zero with an 
average V re i < 0. This implies that the neutron veloc- 
ity is smaller than the fragment velocity (at half target 
thickness) . This would be the case if the neutrons were 
emitted after the 26 O traveled through a portion of the 
target decreasing its velocity. 

It is important to understand the calibration, reso- 
lution, and accuracy of the neutron and fragment ve- 
locities. The neutron timc-of-flight (and therefore ve- 
locity) was calibrated based on the time-of-flight of the 
gamma-rays produced at the target traveling to MoNA. 
From the width of the gamma peak the relative resolu- 
tion (FWHM/centroid) of the neutron velocity is about 
3%. In comparison to the resolution, the accuracy of the 
neutron velocity is determined from the accuracy of the 
time-of-flight measurement and location of the MoNA 
detector. The neutron velocity accuracy was determined 
to be 0.03 cm/ns at beam velocity (11.8 cm/ns). The 
fragment velocity is determined from the track of the 
fragment through the dipole magnet which is measured 
using two Cathode Readout Drift Chambers (CRDCs). 
The accuracy of the fragment velocity is related to the ac- 
curacy of the magnetic field map and measured position 
of the CRDCs. Reasonable variation of these parameters 
showed the fragment velocity to be 0.02 cm/ns at beam 
velocity. While the resolutions will determine the width 
of the V re i distribution, the accuracy of the centroid is 
related to the accuracy of the neutron and fragment ve- 
locity measurements. The detector resolutions were in- 
cluded in the Monte Carlo simulation. 

In order to extract a half-life limit (Ti/ 2 ) of 26 O from 
the V re i distribution, the Monte Carlo simulation was 
modified such that the probability distribution for the 
26 decay based on Tx/i was included. Thus, after the 
one-proton reaction occurred within the target (at a ran- 
dom position) the 26 O was propagated for a time t deter- 
mined from the probability distribution before decaying 
into 24 O + n + n. The resulting V re i distributions from 
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0, 4, and 10 ps is compared to 



the experiment in Fig. |3Ja). The V re i distribution with 
T\/2 =0 is unable to reproduce the experimental data. A 
much better fit is achieved with T\/2 = 4 ps, which shows 
a similar shift in the simulation as the experiment. This 
suggests that 26 O did not decay instantaneously but had 
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FIG. 3. (Color online) (a) Experimental V re i distribution from 
the decay of 26 O compared to the Monte Carlo simulation 
where the 2e O half-life is set as ps, 4 ps, and 10 ps. (b) 
Negative log-likelihood (— ln[L]) as a function of the half-life 
from the unbinned fit of the experimental V re i data set. The 
la, 2a, and 3a confidence interval are indicated. 



an appreciable lifetime. The shape of the decay energy 
spectrum (Fig. [1]) would also be affected by the finite life- 
time of 26 O. The Monte Carlo simulations showed that 
significant changes in the Edecay spectrum would be ob- 
served for T1/2 > 10 ps. 

Due to the low statistics of the experiment, the \ 2 
analysis was observed to be dependent on the binning 
of the data. Therefore, a unbinned maximum likeli- 
hood technique was employed to determine the statis- 
tical significance of the results. This procedure is de- 
scribed in Ref. [32^ and was recently used in the analy- 
sis of 27 - 28 F measured with MoNA [2§ [33|. The neg- 
ative log-likelihood (— ln[L]) is plotted as function of 
T1/2 in Fig. Elb). A minimum in — ln[L] is found at 
4.5 ps. The na confidence intervals are calculated as 
ln[L ma J— ln[L] < n 2 /2. The la, 2a, and 3a confidence 
intervals are shown in Fig. EJb). The statistical signif- 
icance of the results indicate that 26 O has a half-life of 



4 




-0.5 0.5 

V rel (cm/ns) 



FIG. 

23, 



4. Experimental V re i distribution from the 
1 0* — > 22 0+n decay. The Gaussian fit (dashed line) is shown 
to guide the eye. 



about 4.5 ps. 

In addition to the statistical significance, it is impor- 
tant to account for possible systematic uncertainties. As 
previously discussed, the accuracy of the neutron and 
fragment velocities was 0.03 and 0.02 cm/ns, respectively. 
This represents a total systematic uncertainty of 0.05 
cm/ns in the V re i distribution, which corresponds to a 
1.7 ps systematic uncertainty, and indicates a finite half- 
life of 26 at 95% confidence level. The systematic uncer- 
tainty was also estimated through examining the neutron 
decay of the first excited state of 23 0* — > 22 0+n, which 
was also populated during the experiment from the 27 F 
beam. Thus, the 9 Be target, Bp of the dipole, MoNA 
configuration, Sweeper detector settings, and calibrations 
were identical to the 26 O measurement and can provide 
an estimate of any unknown systematic errors. Since the 
decay of 23 O should not have a long-lived component, 
the relative velocity spectrum should not be shifted away 
from V re i = (see Fig. 2]). Following the same half-life 
analysis discussed above, the upper limit on la T\ii was 
3 ps for the 23 O distribution. Therefore, the system- 
atic uncertainty was estimated as 3 ps in comparison to 
1.7 ps determined above. The half-life of 26 O is then 
taken as (stat.) ±3(sys.) ps, which gives T X j 2 > 

at 82% confidence level. A new measurement with im- 
proved statistics would allow for both the statistical and 
systematic uncertainties to be reduced. 

Possible systematic effects related to the selection of 
the 24 O events and the application of the 2n causality 
cuts were also investigated. The selection of the 24 O 
events (shown in Fig. 2 of Ref. [13) based on the time- 
of-flight was varied to examine the dependence of the 
V re i distribution (or half- life). The results showed minor 
variations and the half-life remained within the la limit 
shown in Fig. [3J While the causality cuts allow for the 
removal of the large majority of all false 2n events, it is 
important to verify that shift in the V re i distribution is 
not created by the causality cuts. The causality cuts were 
removed from both the experimental and simulated data 
and the shift in the V re i distribution was maintained. 

The results of the lifetime analysis are compared with 
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FIG. 5. (Color online) Predications of Grigorenko et al. for 
the true two-neutron emission width (half-life) as a function of 
the decay energy. The grey region represents the 3a statistical 
limit on the Ti/2 from the current work and the ground state 
resonance energy limit from Ref. [TtJ • This figure was adapted 
from Ref. 



the prediction of Grigorenko et al. [l6[ in Fig. [5] The 
solid grey region is defined from the 3er statistical limit 
on the T1/2 and the ground state resonance energy limit 
of E < 200 keV [13 for 26 O. As shown the agreement 
with the predictions will depend greatly on improving 
the constraints on the energy and the configuration of 
the ground state. For example, a small ^(s) 2 component 
in the 26 O ground state will require the resonance energy 
to be very small for Ti/ 2 = 4.5 ps. 

In summary, a new technique for measuring the life- 
times of neutron unbound nuclei has been presented and 
applied to the case of 26 O. A shift in the relative velocity 
between the 24 and emitted neutrons was observed and 
shown to be related to the lifetime of 26 O. Detailed Monte 
Carlo simulations, in which the half-life of 26 could be 
varied, were compared to the experimental data. The 
extracted 26 half-life was 4.5^'g (stat.) ± 3(sys.) ps. 
This corresponds to 26 O having a finite lifetime with an 
82% confidence level and suggests the possibility of two- 
neutron radioactivity. It appears that, much like the pro- 
ton dripline, the unique structure of the neutron dripline 
nuclei opens the door for observations of new modes of 
radioactivity. Future experimental work is needed to con- 
firm this observation and provide stringent constraints on 
the properties of the 26 ground state. 
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